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ABSTRACT 

Peculiar velocity measurements are the only tool available in the low-redshift 
Universe for mapping the large-scale distribution of matter and can thus be used 
to constrain cosmology. Using redshifts from the 2M-I-I- redshift compilation, we 
reconstruct the density of galaxies within 200 h~^Mpc, allowing for the first time 
good sampling of important superclusters such as the Shapley Concentration. We 
compare the predicted peculiar velocities from 2M-b-b to Tully-Fisher and SNe 
peculiar velocities. We find a value of = 0.431 ± 0.021, sug¬ 
gesting = 0.401 ± 0.024, in good agreement with other probes. The 

predicted peculiar velocity of the Local Group arising from the 2M-b-b volume 
alone is 540 ± 40 km s“^, towards I = 268° ±4°, b = 38° ± 6 °, only 10° out 
of alignment with the Cosmic Microwave Background dipole. To account for ve¬ 
locity contributions arising from sources outside the 2M-I-I- volume, we fit si¬ 
multaneously for and an external bulk fiow in our analysis. We find that an 
external bulk fiow is preferred at the 5.la level, and the best fit has a velocity of 
159 ± 23 km s“^ towards I = 304° ± 11°, 5 = 6° zb 13°. Finally, the predicted 
bulk fiow of a 50 h~^Mpc Gaussian-weighted volume centred on the Local Group 
is 230 ± 30 km s“^, in the direction I = 293° zb 8 °, 5 = 14° zb 10°, in agreement 
with predictions from ACDM. 


1 INTRODUCTION 

Peculiar velocities, i.e. deviations in the motions of galaxies 
from the Hubble flow, are valuable tools which probe the un¬ 
derlying distribution of dark matter, and are in fact the only 
practical means of doing so on large scales in the low red- 
shift Universe. They can be used to constrain the amplitude 
of matter density fluctuations on a range of scales. As the 
amplitude of such fluctuations are themselves cosmology- 
dependent, the analysis of peculiar velocities provides a di¬ 
rect means of testing cosmological predictions. 

In the current standard cosmological paradigm, the ob¬ 
served structure in the Universe is a result of gravitational 
instabilities which grew from density perturbations in an 
otherwise homogeneous background. This gravitational at¬ 
traction of objects to surrounding structure results in pecu¬ 
liar motion, i.e. motion in addition to that resulting from the 
expansion of the Universe. Assuming only mass continuity 
and standard gravitation in an expanding universe, in the lin¬ 
ear regime where these fluctuations in density are small, i.e. 
S{r) = {p — p)/p < 1, peculiar velocities are proportional 


to gravitational accelerations. This relation, as expressed in 
integral form, is as follows: 



where v(r) is the peculiar velocity held, (i(r) is the mass 
density contrast, Um is the cosmological density parameter, 
and where distances are measured in km s“^ {i.e. r — HR, 
where H is the Hubble parameter and R is the comoving 
distance in Mpc). The growth rate of density perturbations, 
is generally parameterized by UJJi, where 7 = 0.55 
for ACDM (Wang & Steinhardt 1998), which has recently 
been shown to be consistent with observations (Hudson & 
Turnbull 2012). 

As the total matter density contrast cannot be observed, 
however, to make use of Equation (1), an assumption must 
first be made as to how observed galaxies trace the under¬ 
lying total matter. Assuming linear biasing holds on large 
scales, Sg = bS, where b is the linear bias factor, and where 
we have used the subscript “g” when referring to galaxies. 
Rewriting Equation (1) in terms of the observable density 
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contrast of galaxies in the nearby Universe, the proportional¬ 
ity factor between gravitational acceleration and peculiar ve¬ 
locity is then [3 = f /b. Thus, if linear theory holds, by com¬ 
paring measured peculiar velocities to those predicted by 
the distribution of galaxies in redshift surveys, one can con¬ 
strain cosmological parameters through the measurement 
of (3. Furthermore, under the assumption of linear biasing 
cr8,g = has, where as is the root mean square density fluc¬ 
tuations on an 8 h~^Mpc scale. 

By measuring as,g directly from the redshift data, one 
can eliminate b and constrain the degenerate cosmological 
parameter combination fas = l^as^g. Using this method 
of velocity-velocity comparison, a number of recent studies 
have eonstrained this parameter combination (Pike & Hud¬ 
son 2005; Davis et al. 2011; Branchini et al. 2012; Turnbull 
et al. 2012). 

In addition to constraints placed on these parameters, 
when full sky surveys are used one can compute the ve¬ 
locity of the Local Group (hereafter LG) arising from the 
volume under eonsideration as predicted by linear theory. 
While it has long been assumed that the dipole in the Cos¬ 
mic Microwave Background (CMB) temperature map is a 
Doppler effect due to the Sun’s motion, this has only been 
proved recently, via the aberration of the CMB tempera¬ 
ture anisotropies (Planck Collaboration et al. 2013b). The 
motion of the Sun with respect to the Galaxy and of the 
Galaxy with respect to the barycentre of the LG are well 
known (e.g. Courteau & van den Bergh 1999), and when 
combined with the Sun’s motion with respect to the CMB 
(Hinshaw et al. 2009) allows a determination of the motion 
of the LG with respect to the CMB: 622 zb 35 km s“^ in the 
direction I = 272° zb 3°, b — 28° zb 5°. A deviation from 
the predicted value with that derived from the CMB dipole 
would presumably arise from sources beyond the survey and 
would thus have implications for large scale structure. As the 
tidal field falls off as to first order one can model addi¬ 
tional velocity contributions to the LG arising from sources 
outside the survey volume as a dipole, with the magnitude 
of this dipole, or residual bulk flow, itself being a test for 
cosmological models. Past studies constraining cosmologi¬ 
cal parameters through comparison of predicted motion of 
the LG using linear theory and that derived from the CMB 
include Erdogdu et al. (2006), and Bilicki et al. (2011). Re¬ 
cently such an analysis has been extended to the non-linear 
regime using a novel orbit-reconstruction algorithm to pre¬ 
dict motions of nearby objects (Lavaux et al. 2010). 

In addition to reconstructing the motion of the Local 
Group, one can explore the bulk motion of a large volume 
(typically the mean velocity of a 50 h~^Mpc Gaussian- 
weighted window) as such motion probes the amplitude of 
matter power spectrum on large scales. As previous studies 
in this vein have found hints of excess power on large scales 
(cf. Watkins et al. 2009), there has been interest in perform¬ 
ing such analyses using peculiar velocity surveys (e.g. Turn- 
bull et al. 2012; Hong et al. 2014). 

In this work we explore the methods used to self- 
consistently reconstruct the real space density field from red¬ 
shift space while quantifying any biases intrinsic to this re¬ 
construction method. We then use the recently constructed 
large full-sky 2M-1-I- eatalogue composed of 69,160 galaxy 
redshifts to measure /3 to high precision and constrain fas. 
We further explore the growth of the LG’s dipole as pre¬ 


dicted by linear theory arising from structures within 2M-I-+, 
in addition to computing the bulk flow arising from this sur¬ 
vey. 

This paper is organized as follows: in §2 we briefiy re¬ 
view construction of the 2M-I-1- catalogue, we discuss ac¬ 
counting for incompleteness and functional dependence of 
galaxy bias on luminosity when computing the density field, 
and we outline details of the reconstruction procedure. In §3 
we explore the growth of the LG velocity amplitude arising 
from 2M-I-I- as a function of survey depth. In §4 we discuss 
peculiar velocity surveys used, the methods used in com¬ 
paring predicted velocities to measured velocities, as well 
as the results obtained from such analyses. We discuss and 
compare our results to those from recent literature in §5, and 
conclude in §6. 


2 DENSITY FIELD RECONSTRUCTION 

Redshift surveys measure positions of objects in redshift- 
space. As Equation (1) requires real-space positions, we 
must first map observed redshifts to real-space distances. 
The observed redshift, Zobs, is related to cosmological red¬ 
shift resulting from the expansion of the Universe, ^cos, and 
that resulting from peculiar velocities, Zpec, through 

(1 “b -S^obs) — (1 “b ^cos )(l+^pec), (2) 

where for non-relativistic peculiar motions czpec — ^pec- 
Note that in the above if Zobs is corrected to some frame of 
reference, such as the CMB, then Upec is the peculiar ve¬ 
locity with respect to that frame. The comoving distance in 
the low-redshift Universe is then related to the cosmological 
redshift through 

HoR ~ C (^^cos - , (3) 

(Peebles 1993) where R is the comoving distance, Hq is 
the local value of the Hubble parameter, and where is 
the local value of the deceleration parameter, which is given 
by ^ = Dm/2 — Da for a fiat ACDM Universe. As real- 
space positions are dependent on peculiar velocity predic¬ 
tions, which are themselves dependent on real-space posi¬ 
tions, mapping redshifts to real-space must be done with 
care. 

In this section, we discuss the 2M-I-I- redshift catalogue, 
and outline the procedure used in reeonstructing comoving 
positions of the galaxies therein. §2.1 reviews construction 
of 2M-I--I-, and in §2.2 we outline how galaxy weights were 
computed to aceount for the fact that 2M-1-I- is magnitude 
limited. In §2.3 we briefiy review the choice of smooth¬ 
ing kernel. We discuss the procedure used to normalize the 
smoothed density field to the same effective bias in §2.4. In 
§2.5 we outline the iterative scheme used to recover real- 
space positions from redshift-space, and in §2.6 we take a 
cosmographic tour through the recovered density field. 

2.1 2M-H- Redshift Compilation 

The integral in Equation (1) is over all space. In reconstruct¬ 
ing the velocity field of the local Universe, therefore, clearly 
one would like a redshift survey that is very deep and as 
close to all-sky as possible. Two such catalogues which have 
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been used extensively in the past include the sparsely sam¬ 
pled IRAS Point Source Catalogue Redshift Survey (PSCz, 
Saunders et al. 2000), and more recently, the shallower but 
more densely sampled Two-Micron All-Sky Redshift Survey 
(2MRS, Huchra et al. 2012). In this work we use a superset 
of 2MRS, dubbed 2M-I-I-, constructed by Lavaux & Hudson 
(2011). This sample has greater depth than 2MRS, and supe¬ 
rior sampling than PSCz. The photometry is from the Two- 
Micron-All-Sky-Survey (2MASS) Extended Source cata¬ 
logue, (2MASS-XSC, Skrutskie et al. 2006), an all-sky sur¬ 
vey in the J, H and bands. Redshifts in the band 
of the 2MASS Redshift Survey (2MRS) are supplemented 
by those from the Sloan Digital Sky Survey Data Release 
Seven (SDSS-DR7, Abazajian et al. 2009), and the Six- 
Degree-Field Galaxy Redshift Survey Data Release Three 
(6dFGRS, Jones et al. 2009). Data from SDSS was matched 
to that of 2MASS-XSC using the NYU-VAGC catalogue 
(Blanton et al. 2005). As 2M-I-I- draws from multiple sur¬ 
veys, galaxy magnitudes from all sources were first recom¬ 
puted by measuring the apparent magnitude in the band 
within a circular isophote at 20 mag arcsec”^. Following a 
prescription described in Lavaux & Hudson (2011), magni¬ 
tudes were then corrected for Galactic extinction, cosmolog¬ 
ical surface brightness dimming and stellar evolution. After 
corrections the sample was limited to K 2 M++ ^ 11.5 in 
regions not covered by 6dFGRS or SDSS, and limited to 
K 2 M++ ^12.5 elsewhere. 

Other relevant corrections which were made to this cat¬ 
alogue include accounting for incompleteness due to fibre- 
collisions in 6dF and SDSS, as well as treatment of the zone 
of avoidance (ZoA). Incompleteness due to fibre-collisions 
was treated by cloning redshifts of nearby galaxies within 
each survey region as described in Lavaux & Hudson (2011). 

In treating the ZoA, for Galactic longitudes in the range 
[30°, 330°], lower latitudes (\b\ <5°) were first masked and 
then cloned with the redshifts from 2MRS in an equal-area 
strip just above the missing northern (0° < 5 < 5°) Galac¬ 
tic strip, and, for the negative Galactic latitudes in the south, 
a strip below that was cloned. Near the Galactic centre, for 
longitudes in the range [-30°, 30°], the wider Galactic lati¬ 
tude strip \b\ <10° was filled with the redshifts from 6dF- 
GRS in a similar way. A histogram of distances is shown in 
Figure 1. 


2.2 Luminosity Function & Galaxy Weights 

Before using the catalogue to construct the density field, we 
must first account for survey incompleteness. In this section 
we provide a summary of the method used to obtain the lu¬ 
minosity function fit to the catalogue; this luminosity func¬ 
tion is in turn used to in the weighting scheme employed to 
account for incompleteness. For a complete description of 
these calculations as applied to this catalogue see Lavaux & 
Hudson (2011). The luminosity function used to character¬ 
ize the dataset is the Schechter function (1976), which when 
written in terms of absolute magnitudes is given by: 

T>(M) = 

0.41og(10)n*10° ‘‘'^+"^^^*“^^exp ^_ioO'4(Jw*-m) j ^ 

( 4 ) 


where n* is the density normalization, M* is the abso¬ 
lute magnitude break, and a is additional power-law pa¬ 
rameter to be determined. Schechter function parameters are 
computed using likelihood formalism, where the product of 
all conditional probabilities of observing galaxies intrinsic 
magnitudes is maximized given their redshifts, Schechter 
parameter values, and survey completeness at their specified 
angular positions and distances. 

Before computing weights we first discuss the differ¬ 
ent ways by which we can model the galaxy density con¬ 
trast. As we cannot compute the mass density contrast of 
observed galaxies directly, we must use either the number- 
density of galaxies, or their luminosity-density in comput¬ 
ing Sg. In the context of linear biasing, our goal is to create 
a galaxy density contrast field which most closely traces the 
underlying total mass density contrast. Although luminosity- 
density may be a better proxy for stellar mass, and thus for 
the underlying mass distribution of dark matter, we will con¬ 
sider both schemes in this work. 

To account for incompleteness, galaxies are weighted 
according to a common prescription similar to that of Davis 
& Huchra (1982). In the case of a number-density scheme, 
observed galaxies are weighted to account for the number 
of galaxies not observed at a given distance due to the mag¬ 
nitude limit of the survey. When using the galaxy number- 
density for a single homogeneous redshift survey, galaxies 
are weighted by 


Nr \ 

w (r) 


iVa^ 


-Yobserved ) 




(5) 


when 4:7Vdj^fmin > Lmin, and unity otherwise. The flux 
limit, /min, corresponds to a Ks band apparent magnitude 
limit of 11.5 for galaxies drawn from 2MRS and 12.5 other¬ 
wise. The luminosity Lmin used above corresponds to a Ks 
band absolute magnitude of —20. Computed weights used 
in this work additionally account for the inhomogeneous in¬ 
completeness of 2M-I-I- and a complete description of their 
determination can be found in Lavaux & Hudson (2011). 

In this paper, by default, we will use luminosity-density 
to compute the galaxy density contrast. The weight assigned 
to each galaxy’s luminosity is again based on the fraction 
of the total luminosity expected, given the magnitude limit 
of the survey, to the luminosity one expects to observe at 
a given distance. Thus, for a single homogeneous redshift 
survey, galaxy luminosities are weighted by 


L / \ -Laverage 

w [r) — 


^observed ) 


L T>(L) dL 


( 6 ) 


The calculated luminosity weights obtained for the best val¬ 
ues of cosmological parameters are shown in Figure 2. 


2.3 Smoothing 

In order to use Equation (1), the density field must first be 
sufficiently smooth for linear theory to apply. The optimal 
scale on which to smooth the data was determined by com¬ 
paring velocities from an N-body simulation to predictions 
obtained through linear theory using different smoothing 
lengths. Smoothing the density contrast with a 4 h~^Mpc 
Gaussian was found to be the best compromise in minimiz¬ 
ing the scatter in predicted velocities vs. simulation veloci- 
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Figure 1. Left: All galaxies in 2M++ with measured redshifts; blue galaxies are nearest, red are farthest. The 2MRS region {K 2 M++ ^ 11.5) 
with lower density of galaxies is apparent. The Galactic centre is in the centre of the plot, and Galactic longitude increases to the left. Right: 
Histogram of galaxies in catalogue as a function of distance (bin width of 10/i“^Mpc). 



Figure 2. Computed luminosity weights obtained from the proce¬ 
dure outlined in §2.2. As 2MRS is magnitude limited to Ks < 11.5, 
weights rise more sharply than for SDSS and 6dF galaxies which 
are magnitude limited to iCg < 12.5. The scatter in the weights at 
a given distance arise from varying levels of redshift completeness 
across the sky. 


ties, while simultaneously returning an unbiased slope in the 
comparison of observed nonlinear velocities (from the sim¬ 
ulation) with the linear theory predictions from a smoothed, 
reconstructed density field of halos. These comparisons are 
discussed in greater detail in Appendix A. 


2.4 Accounting for Magnitude Dependence of 
Galaxy-Matter Bias 

For magnitude limited surveys such as 2M++, the mean 
luminosity of observed galaxies increases with depth. As 
galaxy-matter bias has been found to increase with luminos¬ 
ity, this means that objects observed at higher redshift are on 
average more biased than those observed nearby. In this sec¬ 
tion we account for this effect by rescaling the density field 
to the same effective bias. We do so using the bias model of 
Westover (2007) in which bias is a function of luminosity. 


By comparing the correlation function of 2MASS volume- 
limited subsamples, converting the binned absolute magni¬ 
tude to a luminosity and defining h/h* — 

Westover (2007) found h/h* = (0.73 ± 0.07) + (0.24 ± 
0.04)L/L*, where h* is the bias of an L* galaxy. This result 
is consistent with that of Norberg et al. (2001) and Tegmark 
et al. (2002), who did similar analyses using projected corre¬ 
lation functions of 2dFGRS and the SDSS power spectrum, 
respectively. 

Since bias is a function of luminosity, this means that 
a naive computation of the density field using a magnitude 
limited survey would lead to a larger effective bias at larger 
distances. As our end goal is to compare predicted velocities 
with measured velocities and determine //5, we must first 
correct the density contrast field by normalizing the field to 
the same effective bias. The effective number-weighted bias 
is computed as follows: 






{r)b*. 


(7) 


Using a luminosity-weighting scheme in computing the den¬ 
sity contrast, the effective luminosity-weighted bias is given 
by: 


bifsir) 


JZ2f^,JiL)L<i>iL)dL 


■ip^{r)b*. (8) 


Using the functional form of 6(L) quoted above from 
Westover (2007), this normalization procedure was applied 
to the 2M-I-1- density contrast fields. Both the number- 
weighted and luminosity-weighted effective bias for a = 
—0.85, M* = —23.25, and a magnitude limit of 12.5 are 
shown in Figure 3. We can then rewrite Equation (1) as 


v(r) 



(9) 


where we have defined = /(Q)/5*, and Sg{r) = 
b*5{r) = 5g{r)/ip{r). 
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Figure 3. Number-weighted (lower) and luminosity-weighted (up¬ 
per) effective bias as a function of distance for the most signifi¬ 
cant la deviation of parameters from the scaling relation 6/6* = 
(0.73 ± 0.07) + (0.24 ± 0.04)L/L*. Plots are obtained using the 
parameter values a = —0.85 and M* = —23.25 for the Schechter 
luminosity function. 

2.5 Reconstruction Procedure 

The weighted galaxies from 2M++ have measured redshifts 
cz, and not precise distances r. However, application of 
(1) requires distances as opposed to redshifts. We refer to 
the inverse problem of determining the positions from red¬ 
shifts as “reconstruction.” Reconstruction was accomplished 
via an iterative procedure modeled on that of Yahil et al. 
(1991). Objects were first grouped using the “Friends-of- 
friends” algorithm (Huchra & Geller 1982), and then placed 
at the mean of their group redshift distance to suppress the 
“Fingers-of-God” effect. Gravity was then “adiabatically” 
turned on by increasing /3* = /(Om)from 0 to 1 in steps 
of 0.01. The reconstruction took place in the LG frame, and 
on each iteration the following steps were taken: 

(i) A Schechter luminosity function is fitted to the data 
using the likelihood formalism discussed in §2.2. The LF 
is in turn used to compute either the luminosity or number 
weights following the procedure discussed in §2.2. 

(ii) Galaxies from 2MRS with distances greater than 125 
h~^Mpc are assigned a weight of zero. Galaxy properties, 
including newly computed weights are then cloned to ac¬ 
count for incompleteness and fill the ZoA as described in 
§ 2 . 1 . 

(iii) Number weighted galaxies or their weighted lumi¬ 
nosities within 200 /i“^Mpc are then placed on a grid. 

(iv) The density constrast field is then computed and nor¬ 
malized to the same bias, 6*, as described in §2.4. The field is 
in turn smoothed with a Gaussian kernel of width 4 h~^Mpc. 

(v) Using Equation (9) the density contrast field is then 
used to obtain predicted peculiar velocities for all objects in 
the catalogue. 

(vi) In conjunction with measured redshifts, predicted pe¬ 
culiar velocities projected on to the line-of-sight are then 
used to predict comoving distances using Equations (2) and 

(3). 

(vii) The previous five predictions for a galaxy’s distance 
are then averaged to suppress oscillations arising from triple 
valued regions, i.e. regions of high-density near which there 


are multiple solutions for distance given redshift (discussed 
further in Appendix B). The averaged distance is in turn as¬ 
signed to the galaxy and used to recompute the galaxy’s ab¬ 
solute magnitude. Computed distances and magnitudes are 
then used in the subsequent iteration. 

Catalogues containing updated distances, luminosities and 
weights were saved at each iteration (corresponding to in¬ 
creasing values of (3) in addition to the computed density 
and velocity fields. 

It should be noted that this iterative reconstruction pro¬ 
cedure was found to be unbiased in determining the best fit 
value of when the full analysis was run on an N-body 
simulation using a ACDM cosmology (see Appendix A for 
more details). 

2.6 Cosmography 

Eigure 4 shows the luminosity-weighted density field of the 
Supergalactic Plane for /3* = 0.43, smoothed with a 4 
/i“^Mpc Gaussian kernel. The incomplete coverage due to 
the lower magnitude limit of 2MRS is clearly visible in this 
figure beyond SGX ~ 125 h~^Mpc. The most prominent 
overdensity in this plane is the Shapley Concentration lo¬ 
cated at (SGX, SGY) - (-125, 75) /i^^Mpc. Other notable 
structures in Eigure 4 include the Virgo Supercluster directly 
above the LG, the Hydra-Centaurus Supercluster at (-40, 
20) /i“^Mpc, and the Perseus-Pisces Supercluster (40, -30) 
/i“^Mpc. Additional slices through SGZ are shown in Eig- 
ure 5, though smoothed on a 7 /i“^Mpc scale to enhance the 
contrast of large overdensities. Eor instance, Horologium- 
Reticulum Supercluster can be readily observed at SGZ ~ 
-112 /i-^Mpc, SGX - -70 /i-^Mpc, SGY - -140 /i^^Mpc. 
The supergalactic plane is also shown with this smoothing 
for comparison. 

3 2M++ PREDICTED PECULIAR VELOCITY OF 
THE LOCAL GROUP 

The velocity of the LG as predicted by linear theory for an 
ideal distance-limited catalogue is given by: 

Q rRmax 

VLG = d^r'S^ir') ^ + Vext , (10) 

where Vext encapsulates contributions from beyond i^max, 
and to first order, can be approximated as a dipole, or “resid¬ 
ual” bulk flow. For a realistic flux-limited catalogue, we do 
not detect a continuous distribution of matter but a finite 
number of galaxies. As a result, our estimate of the veloc¬ 
ity of the LG using 2M-1-I- is subject to shot noise. To es¬ 
timate the effect of shot noise on our predicted motion of 
the LG, we computed the standard deviation in each of the 
components of 500 bootstrap samples. The shot noise in the 
amplitude of the LG’s motion from this analysis was found 
to be 56 km s“^. 

Under the assumption that the observed CMB dipole 
arises from the motion of the LG, there has been much de¬ 
bate as to the structures sourcing this motion. Most recent 
work has made use of 2MASS-XSC or 2MRS in reconstruct¬ 
ing the motion of the LG, and there has yet to be consen¬ 
sus on the distance at which the LG’s motion coincides with 
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Figure 4. The Supergalactic Plane (SGZ= 0) of the 2M++ luminosity-weighted galaxy density contrast held, reconstructed with /3* = 0.43 
smoothed with a Gaussian kernel of radius 4 h“^Mpc. The dashed contour is S* = —0.5, the bold white contour is S* = 0, and successive 
contours thereafter increase from 1 upwards in steps of 3. The Galactic plane runs roughly along the SGY= 0 axis. The Shapley Concentration 
is located at (SGX, SGY) ~ (-125, 75) h“^Mpc, the Virgo Supercluster directly above the LG, the Hydra-Centaurus Supercluster at (-40, 20) 
h“^Mpc, and the Perseus-Pisces Supercluster is at (40, -30) h“^Mpc. The density held is shallower at positive SGX because this region is 
only covered by the 2MRS, whereas the rest of the plane is covered by the deeper 6dFGRS and SDSS. 


that derived from the CMB. Erdogdu et al. (2006) argue that 
more than 70% of the LG’s motion results from structures 
within 50 h~^Mpc, such as Hydra-Centaurus Supercluster. 
While others argue for convergence at distances greater than 
~120 /i“^Mpc such as Lavaux et al. (2010) using Monge- 
Ampere-Kantorovich orbit-reconstruction method, or Bil- 
icki et al. (2011), who explored the convergence of the 
2 MASS dipole moment of the angular distribution of galax¬ 
ies as a function of the limiting flux of the sample. As 2M-I-I- 
is a superset of 2MRS and contains redshift measurements 
up to a magnitude limit of iCg = 12.5, this data is well suited 
to examine the influence of structures beyond 120 h~^Mpc 
on the velocity of the Local Group. 

Using density flelds for different values of /3* that 
were obtained throughout the iteration procedure discussed 
in §2.5, the velocity held of increasingly larger concentric 
spheres centered on the LG was computed. The direction 
and amplitude of the LG velocity as a function of distance 
for different values of was then obtained. The amplitude 


of the LG velocity as predicted by 2M-1-I- for different val¬ 
ues of is shown in Figure 6a. The expected agreement 
between predictions using linear theory for a survey of a cer¬ 
tain depth with values derived from the CMB are plotted for 
comparison. This conditional probability assumes a ACDM 
WMAP9 cosmology. A derivation of this conditional veloc¬ 
ity can be found in Appendix A of Lavaux et al. (2010), 
and is based on Lahav et al. (1990). The convergence of the 
direction of LG velocity with that derived from the CMB 
(I = 272° ±3°, 6 = 28° ±5°) is plotted in Figure 6b. For our 
best fit value of /7* = 0.43 from §4 below, the misalignment 
is 10°. This misalignment angle is signiflcantly better than 
those found by past studies using the shallower 2MRS, such 
as the 21° misalignment found by Erdogdu et al. (2006), 19° 
found by Bilicki et al. (2011), or the ~ 45° misalignment 
found by Lavaux et al. (2010). 
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Figure 5. The 2M-1-1- density field at various slices of SGZ. The density field was smoothed with a 7 h~^Mpc Gaussian kernel. The SGZ 
= —110h~^ Mpc panel shows the Horologium-Reticulum supercluster. The dashed contour is S* = —0.5, the bold white contour is S* = 0, 
and successive contours thereafter increase from 1 upwards in steps of 2. 


4 PECULIAR VELOCITY COMPARISONS 

In principle it is possible to constrain cosmological parame¬ 
ters by comparing the growth of predicted motion of the LG 
with its observed motion. This is difficult in practice because 
one expects contributions from sources beyond the survey 
limit, which we have modeled here as a residual dipole Vext • 


There is therefore a degeneracy between ^*and Vext param¬ 
eters appearing in equation (10). Very early studies of the 
growth of the LG gravity dipole observed the flatness at large 
radii from data similar to that shown in Figure 6, therefore 
assumed that f4xt was negligible, and hence solved for /3*. 
This procedure, however, leads to a value of [3 which is bi¬ 
ased high, since Vext and the integral in (10) are correlated. 
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Figure 6. Left: Growth of LG velocity amplitude as predicted by linear theory for successively larger concentric spheres. The non-monotonic 
curves correspond to the 2M-I-I- predicted growth of the LG amplitude for different values of as indicated in the legend. The solid purple 
curve is for our best fit value /3* = 0.43. The smooth thick grey line corresponds to the expected velocity amplitude for a survey of a depth 
indicated by the x-axis in a ACDM WMAP9 cosmology. The thin grey lines indicate the 68% uncertainties due to cosmic variance. The shaded 
cyan band corresponds to the velocity of the Local Group and its 68% uncertainties as inferred from the CMB dipole. Right: The misalignment 
of the LG predicted direction of the velocity arising from 2M++ with that derived from the CMB (I = 272°, h = 28°). The smooth grey curve 
indicates the 68% uncertainty on the alignment expected in the WMAP9 cosmology. 


An alternative approach is to constrain f4xt by assuming 
a cosmological model of density fluctuations, but this then 
makes the exercise model-dependent. 

A better approach is to break the degeneracy by mea¬ 
suring peculiar velocities of galaxies or groups other than 
the LG. In this section, we discuss comparisons between the 
predicted and observed motions of two samples of galaxies, 
derived from the Tully-Fisher (TF) relation and from Type 
la supernovae. 

4.1 Peculiar Velocity Surveys 

4.1.1 SF1+ + 

SFI-i-i- (Springob et al. 2007) builds primarily on Spiral 
Cluster I-band (SCI) and Spiral Field I-band (SFI) samples 
and uses a mixture of 21 -cm line profile widths and opti¬ 
cal rotation curves in determining the I-band Tully-Fisher 
(TF) relation from a subset of 807 galaxies in the fields of 3 1 
clusters and groups Masters et al. (2006). From the derived 
TF relation they in turn determine the peculiar velocities of 
5780 galaxies. Upon removing galaxies without high-quality 
width measurements and those that are located beyond the 
volume covered by 2M-I-I-, SFI-i-i- can be divided in to two 
subsets of 2583 field galaxies and 735 galaxy groups. 

As noted by Davis et al. (2011), the SFI-i-i- TF relation 
has a kink in the faint end (M > — 20 ), and an asymmet¬ 
ric distribution of outliers about the expected velocity width 
parameter 77 = log(IU) — 2.5 (Springob et al. 2007). As we 
will be fitting for the inverse TF relation, we will account for 
the outliers and deviation from linearity of the relation by 
excluding galaxies with redshift-distance magnitudes fainter 
than -20. We then iteratively compute the TF relation param¬ 
eters and remove those with a velocity width that deviates by 
more than 0.2 in 77 (3.8cr) from the relation, until derived pa¬ 
rameter values converge. Selection on both magnitude and 
velocity width resulted in the rejection of 503 field galax¬ 
ies and 137 galaxy groups. Furthermore, when comparing 


predicted velocities from 2M-I-I- with those from SFI- 1 - 1 -, the 
remaining objects which were found to differ by more than 
3.5cr with all velocity fields obtained through the reconstruc¬ 
tion procedure were rejected (0.6%). The final sample was 
composed of 2067 field galaxies and 595 galaxy groups. The 
typical or characteristic depth of the sample can be quanti¬ 
fied by a weighted mean distance, where the weights are the 
inverse square of the uncertainties. This yields depths of of 
42 h~^Mpc and 25 h~^Mpc for field and group samples re¬ 
spectively. 

For the TF relation, we will perform the fit two dif¬ 
ferent ways. The first is a direct maximum-likelihood fit to 
the observed linewidths: the VELMOD method of Willick 
et al. (1997), as described in Section 4.2.1 below. The second 
method uses the estimated distances as given by Springob 
et al. (2007), but corrected for the fact that their peculiar ve¬ 
locities were obtained under the assumption that czohs = 
HoRF Upec. Specifically, we use the analytic relation from 
Equation (2) to obtain velocities from measured positions. 

4.1.2 First Amendment Supernovae 

The Eirst Amendment (Al) catalogue Type la Supernovae 
(SNe) datasets compiled by Turnbull et al. (2012). Al is 
composed of SNe within 200 /i“^Mpc and draws 34 SNe 
from Jha et al. (2007), 185 from Hicken et al. (2009) and 26 
from Eolatelli et al. (2010). Of these 245 SNe, 237 are within 
the volume spanned by 2M-H+, and have an uncertainty- 
weighted depth of 31 h~^Mpc. 

4.2 Velocity-Velocity Comparisons 

Eor the SEI-i-i- subsets we use the distances as determined in 
Springob et al. (2007) which have not been corrected for 
Malmquist bias. We similarly do not use Malmquist bias 
corrected Al distances. Of the comparison methods dis¬ 
cussed below, VELMOD is unaffected by inhomogeneous 
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Malmquist bias, and accounts for homogeneous Malmquist 
bias in the likelihoods. The Forward Likelihood method dis¬ 
cussed accounts for both homogeneous and inhomogeneous 
Malmquist bias in the likelihoods. The simple compari¬ 
son, however, neither aecounts for homogeneous nor inho¬ 
mogeneous Malmquist bias, and as such, the results yielded 
from this analysis are taken to be biased. 


4.2. J VELMOD 


VELMOD is a rigorous maximum likelihood method first 
proposed and implemented by Willick et al. (1997) and de¬ 
scribed further in Willick & Strauss (1998). It is a velocity- 
velocity comparison method used to fit for the TF relation 
parameters (zero-point, slope and scatter) while simultane¬ 
ously fitting for 13. VELMOD takes as inputs TF parameters, 
an object’s redshift and one of the observables (velocity- 
width or apparent magnitude), and maximizes the proba¬ 
bility of observing one given the other. The strength of 
VELMOD analysis is that it neither assumes a one to one 
mapping from redshift space to real space (accounting for 
errors due to triple-valued regions), nor does it require cali¬ 
bration of the TF relation prior to its implementation. 

Forward VELMOD uses the velocity-width to predict 
a galaxy’s apparent magnitude, whereas the inverse method 
uses the apparent magnitude to predict the veloeity-width 
parameter. The forward method is strongly dependent on se¬ 
lection effects, and thus requires a well-modeled selection 
function. The inverse method, however, is much less sensi¬ 
tive to selection effects due to sample selection’s possible 
weak dependence on velocity-width. As the selection func¬ 
tion of SFI-i-i- is rather difficult to model accurately due to 
its being a compilation of various surveys with a range of 
selection criteria, we will make use of the inverse method in 
our analysis. This analysis assumes a TF relation of the form 
f{M) = -brf{M - ainv), where rj = logio(hK) - 2.5, 
M = m — 51og(dL(r)) is the absolute magnitude, c/l is 
the luminosity distanee, and where b[nv, ^inv and arj are the 
slope, intercept and rms scatter of the inverse relation, re¬ 
spectively. 

The conditional probability of observing a measured 
velocity-width of a galaxy with an apparent magnitude, m, 
and an observed redshift, z, is given by: 


where 


P(r]\m, cz) = 


i,cz) = / 
Jo 


P{r], m, cz) 




( 11 ) 


P(if,m,cz)— / dr P(iri,m\r) P{cz\r) , (12) 


P(? 7 , m\r) (X T>(m — /i(r)) x 


S(m^ T], r) exp ( — 


2^2 


(13) 


P[cz\r) — 




: exp 


[cz CZpred] 
2cr2 


(14) 


(1 + ^pred) = (1 + ^cos(r))(l + ^*u(r)/c), (15) 


where Zcos(^) is related to the eomoving distance r = HqR 
through Equation (3), S{m,r],r) is the selection function. 


u is the radial predicted velocity scaled to = 1, is 
the scatter in the aetual velocity compared to the linear the¬ 
ory prediction, and where /i(r) = 5 log r is the distance 
modulus. We take ay = 150 km s“^ based on the tests dis¬ 
cussed in Appendix A. We can then compute the product of 
the conditional probability P{r]\m, cz) over all galaxies and 
minimize the quantity 

Civ = P{r]i\mi, czi) (16) 

i 

for the parameters ^5*, ainv, 5inv, and the three eomponents 
of Vext. 


4.2.2 Forward Likelihood 

We use a maximum likelihood method first described in Pike 
& Hudson (2005) which was developed to compare peculiar 
velocities obtained through SNe surveys while aecounting 
for triple-valued regions. In addition to constraining j3* and 
the three ^mponents of Vext this method can be used to 
constrain h, a nusiance parameter which permits a rescal¬ 
ing of published distances. The forward likelihood method 
maximizes the probability of a galaxy having its observed 
redshift 

noo 

P{cz) = / drP{cz\r)P{r) , (17) 

Jo 

where 

(1 + Zpred) = (1 + Zcos{r,h))(l + I3*u(r)/c), (19) 

where Zcos is given through the relation 

hr = c2cos [^1 - , (20) 

and where 

P(r) ocexp [l+,5|(r)], (21) 

d is the distance as determined by the peculiar velocity sur¬ 
vey, and where ad is the uneertainty is the measured dis¬ 
tance. The product of P{cz) for all objects is then computed, 
from which the quantity £fl = — 2 In P(czi) is mini¬ 
mized. 


4.2.3 Minimization 


In addition to the comparison method diseussed above we 
also perform a simple minimization procedure to deter¬ 
mine the best value of and Vext- For this minimization 
procedure we compare the observed redshift of the object 
with the sum of its measured distance and predicted peculiar 
velocity at that distance, i.e. 


X^(/3*) = LI 

i 


{cZi CZpred) 




2 

V 


( 22 ) 


where Zpred is given by Equation (19). Note that this ex¬ 
pression does not aecount for the effeets of density inhomo¬ 
geneities along the line of sight. The recovered value of is 
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affected by inhomogeneous Malmquist bias and results from 
this method are thus expected to be biased high as a result. 
Nevertheless, the statistic is useful to assess goodness- 
of-fit and so is included here. 


4.3 Results 

Key results obtained through velocity-velocity comparison 
methods discussed in §4.2 are summarized in Table 1. To de¬ 
termine the value of /3* at which £fl and arc minimized, 
a cubic function was fit to resultant data. Comparison of the 
TF relation constants (ainv, ^inv) obtained to those found by 
Masters et al. (2006) are shown in Table 2. The best fit value 
for all parameters, including and its errors were obtained 
from 500 bootstrap samples of both 2M-I-I- and the peculiar 
velocity datasets. Through bootstrap analysis it was found 
that peculiar velocity datasets and 2M-I-I- contributed essen¬ 
tially equal amounts to overall parameter errors. The value 
obtained for the best fit residual bulk fiow, Vext, is given 
in Table 3 along with the bulk fiow (50 h~^Mpc Gaussian- 
weighted mean of the velocity field) and the predicted veloc¬ 
ity of the LG arising from 2M-I-I-. A comparison of with 
and without the residual bulk fiow results in a difference of 
34 for the 3 degrees of freedom, the residual bulk flow model 
is thus preferred at the 5.1 cr level. 

For a qualitative illustration of the agreement of pre¬ 
dicted velocities arising from 2M-I-I- with the published val¬ 
ues from SFI-i-i- and Al we have plotted the projected line- 
of-sight (LOS) velocity within a 30° cone centered on the 
Shapley Supercluster and Hydra Supercluster in Figure 7. 
These two structures were chosen as they lie in approxi¬ 
mately the same direction as Vext • From this figure it is ap¬ 
parent that predictions do in fact follow the trends observed 
in measured velocities. Furthermore, it is apparent that addi¬ 
tion of Vext to predicted velocities seems to provide better 
agreement, suggesting that such a residual bulk fiow is in 
fact warranted. 

In combination with a measurement of crj^g, 13* can 
be used to constrain the cosmology dependent (and sur¬ 
vey independent) degenerate parameter combination fag = 
/3*cr8,g. To measure crj^g from 2M-I-I- we follow a similar 
prescription to that of Efstathiou et al. (1990), we compute 
c^8,g = (^8,g(^)/'0(^)) using counts in cells within radial 
shells. Using this maximum likelihood scheme we obtain the 
value cr8,g = 0.99 zb 0.04, where the errors quoted are de¬ 
rived from the scatter among different shells, as this value 
was found to be more conservative then the formal errors 
from the likelihood analysis. This value is in good agreement 
with that found by Westover (2007) of = 0.98 zb 0.07 
obtained by fitting projected correlation functions to 2MRS 
galaxies within the magnitude range containing L* galaxies, 
Le. —23.5 < Ks < —23.0. The product of the growth factor 
and non-linear ag is thus fag = 0.427 zb 0.026. By adopt¬ 
ing the value of Um = 0.3, we can transform our non-linear 
value of ag to a linearized value following the prescription 
of Juszkiewicz et al. (2010). We in turn obtain the constraint 
/cr8,iin = 0.401 zb 0.024. It is important to note that lin¬ 
earization is only weakly dependent on the adopted value of 
Qm (flm = 0.266 results in fag^un = 0.398 zb 0.024). We 
compare these results with those obtained using independent 
methods in the following section. 


Table 1. Summary of best fit values of /3* using different weighting 
schemes, methods of analysis and peculiar velocity datasets. Re¬ 
sults obtained using luminosity weighting are indicated by (LW), 
whereas those obtained using number weighting are indicated by 
(NW). Unless explicitly indicated, all datasets were used for the 
method mentioned with the exception of Inverse VELMOD which 
used all individual galaxies from SFI-i-i-. 



/3* 

xV(d.o.f.) 

Forward Likelihood (LW) 
Al 

0.440 ± 0.023 


SFI++ Galaxy Groups 

0.429 ± 0.022 

- 

SFI+-I- Field Galaxies 

0.423 ± 0.045 

- 

All 

0.431 ±0.021 

- 

Forward Likelihood (NW) 

0.439 ± 0.020 

- 

Inverse VELMOD (LW) 

0.387 ± 0.048 

- 

(LW) 

0.444 ± 0.026 

2194/2899 

X^ (NW) 

0.442 ± 0.028 

2200/2899 


Table 2. TF relation constants obtained through Inverse VELMOD 
analysis of SFI-i-i- galaxies. Results listed are those obtained using a 
luminosity-weighting (LW) reconstruction scheme. 



U-inv 

^inv 

Masters et al. (2006) 

-20.881 

-8.435 

This Study 

-20.918 ±0.012 

-8.19 ±0.06 


5 DISCUSSION 

5.1 Potential Systematic Effects 

In this section we discuss the various effects which could 
bias our measurement of /S* and our prediction of the motion 
of the LG. Possible contributing factors include: 

(i) Choice of smoothing length used for Gaussian kernel 
may skew our estimate of /S*. 

(ii) Non-linear contributions to velocities from nearby 
small scale structure may systematically bias results. 

(iii) Triple valued regions may result in incorrectly recon¬ 
structed galaxy positions skewing velocity predictions de¬ 
rived therefrom. 

(iv) As the reconstruction is being done in the LG frame, 
at each iteration we are subtracting the motion of the LG 
from all galaxies, our procedure may thus be susceptible to 
the Kaiser rocket effect. 

(v) Sparseness of survey may result in an under¬ 
representation of structures within the survey volume. 

(vi) Structures which lie in the ZoA are not included in 
2M-I-1- and could contribute to the direction and amplitude 
of the motion of the LG as well as influencing predictions of 
peculiar velocities of galaxies near the galactic plane. 

(vii) A non-linear relation between mass and luminosity 
as well as scatter in the underlying relation may infiuence 
predictions. 

We have addressed some of these concerns above, but will 
review them again here for completeness. 

In addressing (i), the width of the Gaussian kernel was 
chosen to be 4 /i“^Mpc as this length was shown to be least 
biased when velocity predictions were compared to those de¬ 
rived from simulations. Furthermore, this smoothing length 
was found to produce minimal scatter in the derived relation. 
We have accounted for (ii) and by performing the full recon¬ 
struction and analysis on data derived from N-body simu- 
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Table 3. The bulk flow and motion of the LG arising from 2M++ for our best value of = 0.43. The best residual bulk flow, Vext, that 
was fitted simultaneously with /3* using the Forward Likelihood is also shown below. The bulk flow was computed by taking a 50 h~^Mpc 
Gaussian-weighted mean of the velocity held corresponding to /3* =0.43. 



(km s ^ ) 

Vy (km s ^) 

Vz (km s ^) 

\\v\\ (km s ^ ) 

longitude (°) 

latitude (°) 

BF 2 M+ + 

-3 ±8 

-72 ± 11 

38± 11 

81 ± 11 

268 ±6 

28± 10 

LG 2 M+ + 

-18 ±27 

-422 ± 41 

328 ± 37 

535 ± 40 

268 ±4 

38 ±6 

Vext 

89 ±21 

-131 ±23 

17 ±26 

159 ± 23 

304 ± 11 

6 ± 13 

BF 2 M++ + Pext 

86 ±22 

-203 ± 26 

55 ±28 

227 ± 25 

293 ±8 

14 ± 10 

LG 2 M++ + Vext 

71 ±34 

-553 ± 47 

345 ± 46 

656 ± 47 

277 ±4 

32 ±6 




Figure 7. Measured velocities of objects from A1 and SFI-i-i- ly¬ 
ing within a 30° cone of Shapley/Great Attractor (upper), and Hy¬ 
dra Superclusters (lower). Velocities are projected on to the line-of- 
sight. Predictions arising from 2M++ for our best value of /3* = 
0.43 for these objects are shown above with and without the pre¬ 
ferred residual bulk flow as connected blue dashed and red solid 
curves, respectively. Values plotted are the error weighted mean of 
velocities within bins of 10 /i“^Mpc. To the left of the origin the 
negative value of the peculiar velocity is plotted so that a bulk flow 
would appear as a constant offset. 


lations. It was found that a value of 150 km should be 
used for the scatter around predictions to account for these 
effects. As for (iii), we used maximum likelihood methods 
which integrate likelihoods along the line-of-sight. As these 
methods account for uncertainty in a galaxy’s position, re¬ 
sults should not be sensitive to a misallocation of a galaxy 
in the event that it is lies within a triple-valued region. The 
value of /3* derived from this analysis was found to be un¬ 
biased when performed on simulations. A more complete 


discussion of the quantification of these systematics through 
analysis of N-body simulations can be found in Appendix A. 

As seen in Figure 6 although the direction of the LG is 
not very susceptible to the value of /3*, the overall amplitude 
varies by ~1200 km s“^ between ^*=0 and ^5*= 0.8. As 
we are doing the reconstruction in the LG frame, an error 
in the estimate of the motion of the LG may result in spuri¬ 
ous distance estimates of objects along this line of motion. 
Putting objects at the incorrect distance may in turn result 
in incorrect object weights, and in turn, incorrect velocity 
predictions. This phenomenon has come to be known as the 
Kaiser “Rocket Effect”, as it was first discussed in Kaiser 
(1987). To account for this effect Strauss et al. (1992) ex¬ 
plored a “Kaiser Fix” to the IRAS 1.2 Jy sample. This “fix” 
amounts to altering the predicted distances of objects to 

r = cz - f • (V(r) - V(0)exp(-rVrK) ^^ 3 ) 
- Vcmb[1 - exp(-r^/r|)]), 

where Vcmb is the velocity of the LG as inferred from the 
CMB dipole, and where vk is 1,000 km s“^ as determined 
by the observed velocity correlation function (Bertschinger 
et al. 1990). Note that this fix assumes that galaxies more 
distant than tk are in fact at rest in the CMB frame, which 
may not be the case; indeed our data suggests otherwise. 
Nevertheless, to estimate the sensitivity of our results to this 
effect, we have implemented the Kaiser fix. We find that the 
final estimate of differs by only 3%, which is small com¬ 
pared to the random errors. 

The impact of survey sparseness on the methodology 
applied in this work has been estimated in the past by Pike 
& Hudson (2005). For the 2MASS catalogue they found 
that under-sampling by 50% produced negligible results on 
their final estimate of /3* (2-3%). Furthermore, we have ac¬ 
counted for the effects of sparse sampling by obtaining our 
quoted results and errors from bootstrap resampling of both 
2M-1-I- and peculiar velocity datasets. In addressing (vi) we 
measured /3* at high latitudes (b > 50°) and found no devia¬ 
tion from previous results beyond that of the random errors. 

When the density field was normalized to the same 
bias, number-weighting and luminosity-weighting schemes 
yielded consistent results. As the was found to be smaller 
for the luminosity weighted result than that obtained with 
number weighting, we hereafter will be quoting the best 
value of as that obtained from the luminosity weighting 
scheme. Results obtained in this paper account for neither 
a non-linear relation between mass and luminosity, nor do 
they account for scatter in mass-luminosity relation. More 
sophisticated models have been proposed, such as the halo- 
model of Marinoni & Hudson (2002), and recently an iter¬ 
ative prescription to reconstruct the density field from the 
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distribution of halos (Wang et al. 2009). We will consider 
implementation of such methods in a future paper. 


5.2 Comparison With Other Results 

As the bias depends on luminosity and possibly morphology 
(c/ Pike & Hudson 2005), the value of 13 obtained is survey 
dependent. As 2M++ draws primarily from 2MASS, how¬ 
ever, a loose comparison can be made with other values of 
/3 obtained therefrom. Most recently, in comparing the clus¬ 
tering dipole of galaxies from 2MASS-XSC to predictions 
from linear theory assuming a ACDM cosmology and con¬ 
vergence of the LG dipole with that derived from the CMB, 
Bilicki et al. (201 1) found (3 — 0.38 zb 0.04. After construct¬ 
ing /3 dependent predictions of peculiar velocities, Branchini 
et al. (2012) in turn estimate {3 by minimizing the scatter of 
predicted 2MASS absolute magnitudes about a universal lu¬ 
minosity function and find {3 — 0.323 zb 0.083. Davis et al. 
(2011) expand the velocity field in spherical harmonics and 
fit the inverse TF relation to SFI-i-i- finding /3 = 0.33zb0.04. 
The TF data span a range of distances, and because the effec¬ 
tive bias changes with distance it is difficult to compare this 
directly with our number-weighted result (3* — 0.44 zb 0.02. 
The characteristic, or error-weighted, depth of all individ¬ 
ual SFI-i-i- galaxies is ~ 32/i“^Mpc. At this distance, the 
typical relative bias of a number-weighted sample with the 
magnitude limit of 2MRS is = 0.93, 

and so applying this to the Davis et al. (2011) result yields 
= 0.31 zb 0.04 

Our value of /erg, un = 0.40 zb 0.02 is in good agree¬ 
ment with those obtained using the same methodology, such 
as that of Turnbull et al. (2012) (0.40 zb 0.07), that of 
Pike & Hudson (2005) (0.44 zb 0.06), as well as with the 
weighted IRAS average of multiple studies reported therein 
(0.40 zb 0.03). As noted above, it is, however, in slight ten¬ 
sion with that found by Davis et al. (2011) (0.31 zb 0.04). 
Note that where necessary the values quoted here have been 
linearized following the procedure discussed above in §4.3. 


Cluster Abundance 


Peculiar Velocities 


0.7 


1 1 


Figure 8. Comparison of fos^un measured results. Values plot¬ 
ted above derived from weak-lensing (Kilbinger et al. 2013, 
CFHTLenS) and cluster abundances [Reichardt et al. (2013, SPT), 
Rozo et al. (2010, MaxBCG), Planck Collaboration et al. (2013a, 
Planck-SZ), Vikhlinin et al. (2009, X-rays)] have assumed a value 
of = 0.3 in mapping constraints to Results obtained 

through previous analyses of measured peculiar velocities are also 
shown [Turnbull et al. (2012, THF), Davis et al. (2011, DNM)], as 
well as from redshift space distortions (Hudson & Turnbull 2012, 
HT). CMB results are from WMAP9 and the Planck Collaboration 
(2015). The horizontal line is the error-weighted mean of all values 
(/fjg = 0.400 ± 0.005), shown here for reference. 


and are at different redshifts, and so their dependence on 
Qm differs in the exponent. To make a quantitative compar¬ 
isons between different results, we adopt Qm = 0.3 (see 
Figure 8). There is some tension between some results e.g. 
Kilbinger et al. (2013) and Planck-SZ (Planck Collaboration 
et al. 2013a) versus Planck CMB temperature (Planck Col¬ 
laboration et al. 2015). The peculiar velocity result presented 
here is consistent with all of these values. 


5.3 Cosmological Implications 

5.3.1 The value of fas,lin 

We can also compare our value fag^nn to constraints placed 
on a degenerate combination of Qm and erg through indepen¬ 
dent means. In particular our value is in excellent agreement 
with a different peculiar velocity probe, namely measure¬ 
ments of f{z)a8{z) at different redshifts via redshift space 
distortions, which yield a best-fit value of fag = 0.40zb0.02 
(Hudson & Turnbull 2012). 

An analysis of second and third-order weak-lensing 
aperture-mass moments measured by CFHTLenS yields 
cr8(fim/0.27)“ ® = 0.79 ± 0.03 (Kilbinger et al. 2013). 
Constraints can also be obtained from the number counts 
and mass of galaxy clusters as measured through X-ray sur¬ 
face brightness (Vikhlinin et al. 2009) and measurements 
of the Sunyaev-Zeldovich (SZ) effect (Planck Collaboration 
et al. 2013a, Reichardt et al. 2013). Finally, we can obtain 
a value for fag from CMB temperature anisotropy from 
Planck = 0.427 ± 0.010; Planck Collaboration 

et al. 2015) and WMAP9 = 0.407 ± 0.029; Hin- 

shaw et al. 2013). The measurements use different methods 


5.3.2 The motion of the LG 

For our best value of we can compare the predicted 
growth of the LG velocity amplitude with the result that 
one would expect to measure using linear theory for a 
ACDM cosmology (conditional on Vcmb). For /3*= 0.431 
zh 0.021, as determined from peculiar velocity comparisons, 
we obtain the prediction for the motion of the Local Group 
arising from 2M-I-I- to be 535 zb 40 km s“^ in the direc¬ 
tion I = 268° zb 4°, 5 = 38° zb 6°, only 10° out of align¬ 
ment with the direction of the motion as inferred from the 
CMB dipole. The residual LG motion is therefore 100 zb 
37 km s“^ in the direction I = 303° ± 36°, b = 34° ± 36°. 
This value is in reasonable agreement with the best fit resid¬ 
ual bulk flow obtained through peculiar velocity compar¬ 
isons in the CMB frame of 159 zb 23 km s“^ in the direc¬ 
tion I = 304° zbll°,6 = 6°zbl3°. Inclusion of this resid¬ 
ual bulk flow with the predicted motion of the LG arising 
from 2M-I-I- results in a total predicted motion of 656 zb 
47 km s“^ in the direction I = 277° zb 4°, 6 = 32° zb 6°, in 
even better agreement with both the amplitude and direction 
of the motion as inferred from the temperature dipole of the 
CMB. 


© 0000 RAS, MNRAS 000, 000-000 













Constraints from 2M++ 13 



Figure 9. Volume-weighted mean of predicted velocity field for 
Gaussian window of increasing scale centered on the Local Group. 
The inferred values from 2M++ with and without the residual bulk 
fiow are shown by the dashed grey line with 68% uncertainties 
in blue hatch, and a dashed black line, with uncertainties in solid 
green, respectively. The predicted root-mean-square velocity for a 
ACDM WMAP9 cosmology is shown as the red solid line, the cos¬ 
mic scatter in the velocity amplitude distribution are shown as red 
dot-dash lines. Bulk fiows in Gaussian-weighted spheres of radius 
40 /i“^Mpc and 50 /i“^Mpc are shown for the results of Hong 
et al. (2014, 2MTF), Turnbull et al. (2012, THF) and Watkins et al. 
(2009, WFH) . The LG motion is also shown, plotted at a radius of 
3 h-^Mpc. 

5.3.3 The residual bulk flow 

We find that the amplitudes and directions of Ve^t fit to each 
of the SFI-i-i- and Al SNe datasets separately are consistent 
with one another. Furthermore, comparing Al with PSCz (of 
comparable depth to 2M-I-+), Turnbull et al. (2012) found a 
residual flow of 14 = 144 ± 44 km s-\ I4 = -38 ± 51 
km s“^, 14 = 20 zb 35 km s“^, in reasonable agreement 
with the values found here of 14 = 89 zb 21 km s“^, I 4 = 
— 131zb23 km s“^, 14 = 17zb26 km s“^. This suggests that 
the residual bulk flow is not an artifact of either the analysis 
or redshift-catalogue and is sourced by structures outside the 
2M-I-I- and PSCz volumes. 

We can also use the 2M-I-I- density field to predict the 
BF and compare this to the BF expected in a ACDM uni¬ 
verse in Figure 9. We have plotted this comparison for the 
Gaussian-weighted mean of the 2M-I-I- velocity field. It is 
apparent from this figure that the resulting bulk flow from 
our analysis is in agreement with that expected for a ACDM 
universe. Combining the cosmic variance in quadrature with 
observational errors, comparison of the measured bulk flow 
of a 100 /i“^Mpc Gaussian with predictions from ACDM 
yield a of 1 -4 for 3 degrees of freedom; clearly the mea¬ 
sured value agrees well with the predicted value from the 
standard cosmological model. 

5.3.4 A large-scale underdensity? 

There have been recent claims that the Local Universe 
(~ 150-200 /i“^Mpc) is under-dense (Whitbourn & Shanks 
2013, Keenan et al. 2013). Such a phenomenon might ac¬ 
count for the discrepancy between the larger value for the 


Hubble parameter when measured locally {z ^ 0) and that 
obtained from studies of the CMB temperature anisotropies. 

Although the majority of 2M-I-I- lies within the sug¬ 
gested underdensity, we have nonetheless explored the 
possibility of a under-dense volume within 2M-I-I-. The 
luminosity-weighted density contrast of 2M-I-I- in shells is 
shown in Figure 10. We have not observed any global sys¬ 
tematic rise in density towards the periphery of the survey. 

To compare our results with others in more detail, note 
that Whitbourn & Shanks (2013) use redshift data from 
three large regions: 6dF-SGC, 6dF-NGC & SDSS-NGC. 
Within z < 0.05, they quote mean density contrast of 
4 = -0.40 zb 0.05, 0.04 zb 0.10 and -0.14 zb 0.05, re¬ 
spectively. For the same z < 0.05 volumes, we find den¬ 
sity contrasts of 4 = —0.17, 0.01 and 0.03 respectively, 
where the density is normalized with respect to the mean 
density within 200 h~^Mpc (z ~ 0.067). Boehringer et al. 
(2014) studied the large-scale densities of X-ray clusters. For 
the 6dF-SGC and 6dF-NGC regions, they find mean cluster 
density contrast of 4i = —0.55 zb 0.10 and 0.02 zb 0.17 
within z < 0.05. However, as they point out, galaxy clus¬ 
ters are highly biased (bd ~ 2.7) and so the corresponding 
mean matter density contrasts are = —0.20 zb 0.04 and 
0.01 zb 0.06. These latter numbers are in good agreement 
with our nearly-unbiased galaxy luminosity results. We con¬ 
clude that, while the 6dF-SGC region may be mildly under- 
dense within < 0.05, there is no evidence for a large-scale 
void. 


5.3.5 Prospects for the future 

There are several upcoming peculiar velocity surveys which 
should dramatically improve the constraints on both and 
fag. Among these is the survey dubbed “Transforming As¬ 
tronomical Imaging surveys through Polychromatic Anal¬ 
ysis of Nebulae” (TAIPAN). Using the UK Schmidt tele¬ 
scope, it is estimated that TAIPAN will acquire ^45,000 
Fundamental-Plane velocity measurements out to a redshift 
of 0.2 (Koda et al. 2013). The next generation of Tully- 
Fisher (TF) peculiar velocity surveys include the Wide- 
field ASKAP L-band Legacy All-sky Blind surveY (WAL¬ 
LABY, Koribalski & Staveley-Smith 2009), and the West- 
erbork Northern Sky HI Survey (WNSHS)^ An HI sur- 
vey acquired using the Australian Square Kilometer Ar¬ 
ray Pathfinder (ASKAP), WALLABY is planned to cover 
37r steradian of sky. Its Northern Hemisphere counterpart, 
WNSHS, is planned to cover remaining tt steradian of the 
sky using the Westerbork Synthesis Radio. It is estimated 
that these surveys will obtain a total of ~32,000 velocity 
measurements, and along with TAIPAN will not only en¬ 
able /c-dependent measurements of fag but will improve 
constraints on this parameter combination at low-redshift 
(z ^ 0.05) to within 3% (Koda et al. 2013). Clearly con¬ 
straints on cosmology through peculiar velocities has a very 
promising future. 


^ http ://www.astron.nl/~j ozsa/wnshs/ 
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Figure 10. The mean density contrast of concentric spherical shells 
of equal volume centered on the Local Group. 


6 CONCLUSION 

Under the assumption of gravitational instability, we recon¬ 
structed the galaxy density field within 200 h~^Mpc using 
2M-I-I-, a composite all-sky redshift redshift survey with high 
completeness. We compared the predicted peculiar veloc¬ 
ities arising from this density field to the measured pecu¬ 
liar velocities of SFI-i-i- and the First Amendment super¬ 
novae dataset using various comparison methods in order 
to measure ^5*= /b*. Using the VELMOD maximum 

likelihood comparison method with SFI-i-i- spiral galaxies 
we simultaneously solved for the inverse TF relation zero- 
point and slope, /3*, as well as a residual bulk flow due to 
sources external to the 2M-I-I- volume, denoted We 
similarly compared peculiar velocities datasets to predic¬ 
tions from 2M-I-I- using a forward-likelihood method in or¬ 
der to constrain (3* and Uext- All methods and data subsets 
yielded consistent values of /3*, with our final result being 
/3*= /(Qm)/^* =0.431 di 0.021. Combining our value of 
/3* with = 0.99 zb 0.04 as measured from 2M-I-I- for 
= 0.43, we in turn measured the parameter combination 
f erg,\in = 0.401 di 0.024. This value was found to be con¬ 
sistent with the majority of results obtained by independent 
means, including those of WMAP9 and Planck. 

For our measured value of we computed the veloc¬ 
ity of the Local Group as predicted by linear theory arising 
from the reconstructed density field. Our value for the veloc¬ 
ity was found to be consistent with the theoretical value that 
would be measured for a survey of this depth in a ACDM 
Universe. Combining our predicted value for the motion of 
LG arising from 2M-I-1- for our best value of with the 
value of Uext obtained through comparing predicted veloci¬ 
ties with peculiar-velocity surveys, we predict a motion of 
the LG to be 660 zb 50 km s“^, towards I = 277° zb 4°, 
5 = 32°zb6°, only 5° out of alignment with the direction 
as inferred from the CMB dipole. Similarly, with addition of 
this residual bulk flow to the 50 h~^Mpc Gaussian-weighted 
mean of the velocity field, we obtained a predicted bulk-flow 
of230zb30km s“^ towards / = 293° zb 8°, 5= 14° zb 10°, 
an amplitude that is consistent with that expected for a 
ACDM Universe. We note, however, that although we find 
the inclusion of Uext is preferred at the 5.la level, it is un¬ 
clear whether this residual bulk flow merely compensates for 


imperfect mapping between luminosity and mass towards 
the periphery of 2M-I-1-, or whether it is due to structures 
outside the 2M-I-I- volume. Future work using more sophisti¬ 
cated biasing schemes may help in answering this question. 

The resulting 2M-1-I- density and peculiar velocity fields 
obtained from this analysis are made available at cos- 
micflows.uwaterloo.ca and cosmicflows.iap.fr. 
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APPENDIX A: TESTS OF THE RECONSTRUCTION 
WITH N-BODY SIMULATIONS 

In this Appendix, we focus on two possible sources of sys¬ 
tematic bias in the reconstructed density and velocity fields. 
First, when constructing the density field from a set of point 
density tracers, it is necessary to smooth to obtain a continu¬ 
ous density field. If the smoothing length is too short, density 
contrasts are high and linear theory is no longer applicable. 
If the smoothing length is too long, then the density con¬ 


trast is suppressed and velocities are underpredicted. Sec¬ 
ond, in “reconstructing” real-space positions from redshift- 
space, the iterative technique discussed in §3 may also in¬ 
troduce a systematic bias in the recovered density field and 
hence in the fitted value of (3. 

Al Effect of smoothing 

Berlind et al. (2000) used N-body simulations to show that 
when predicted velocities derived from smooth density fields 
are compared to measured (unsmoothed) velocities, the re¬ 
covered value of [3 depends on the smoothing. For Gaussian 
smoothing, unbiased results were obtained for a smoothing 
radius between 4 and 5 h~^Mpc. It is interesting to confirm 
these results. 

Here we use an N-body simulation of 512^ particles 
in a 500 h~^Mpc periodic box using GADGET-2 (Springel 
2005). The cosmological parameters of this simulation are as 
follows: Urn = 0.266 , Ua = 0.734, h = 0.71, and where 
each particle is 6.83 x 10^° h~^MQ. From the particle po¬ 
sitions and velocities of the simulation a halo catalogue was 
formed using ROCKSTAR (Behroozi et al. 2013) consisting 
of 693948 halos between 5.5 x 10^^ and 2.2 x 10^^ h~^ 
Mq (between 8 and 31809 particles). Either halo or particle 
positions were smoothed to create a smooth halo or particle 
density field. 

In Fig. Al we plot both the slope of regression between 
smoothed predicted and N-body velocities, as well as the 
scatter about this relation, both as a function of smoothing 
scale. 

For particle velocities compared with linear theory pre¬ 
dictions from the smoothed particle density field, the slope 
is unbiased at a smoothing scale of ~5 /i“^Mpc. The scat¬ 
ter is ~250 km s“^. The results are similar for halos, but 
the scatter is significantly lower: ~150 km s“^. This is be¬ 
cause the particle velocity field is a sum of the motion of the 
halos themselves plus the internal motion of particles with 
respect to the halos. In both cases, the scatter is minimum 
for a smoothing scale of ^4 /i“^Mpc. 

A2 Reconstructing the Halo Distances from Redshifts 

With actual galaxy data, the true distances are unknown but 
redshifts are available. In this paper, we use an iterative re¬ 
construction method to map a galaxy’s position from red- 
shift space to real space. We will refer to these reconstructed 
coordinates as “recon-space.” In the N-body simulation, we 
emulate this by placing halos in redshift-space, and then it¬ 
eratively reconstruct the density field, slowly increasing (3. 

Fig. A2 shows the difference between real-space posi¬ 
tions and recon-space positions. The left panel illustrates the 
absolute displacement in units of km s“^. From this panel, 
it is apparent that in recon-space the majority of halos do in 
fact return to a location close to their real-space positions, 
with an error that is typically ~ 300 km s“^. This is smaller 
than the smoothing scale of ~ 4/i“^ Mpc. However, in prac¬ 
tice this reconstruction error will act as an additional source 
of smoothing. We will return to this point below. 

The right panel shows the displacement of each halo 
from its real-space position normalized by the difference 
between the real-space and redshift-space positions (i.e. by 
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Figure Al. The fitted slope of the regression between N-body observed peculiar velocities and the linear theory predictions as a function of 
smoothing scale is shown in the left panel. In the right panel the scatter about the same regression is plotted. Either particles or halos can be 
used as tracers of the density field, or as tracers of the velocity field. The blue circles and red squares represent halo tracers and particle tracers, 
respectively. The black curve with downwards pointing ‘Y’ show the resulting scatter when the known halo velocities are compared to the 
predictions from the halo density field. The green curve with ‘Y’-shaped symbols show the resulting scatter when the known halo velocities are 
compared to the predictions from the reconstructed halo density field. Note that the reconstruction process shows very nearly unbiased results 
when the field is smoothed with a Gaussian kernel that is 4 /i“^Mpc in length. 



Displacement from simulated "real" position (km s ‘ ) Normalized displacement from simulated "real" position 


Figure A2. Two sample histograms showing the relative displacement of halos in recon-space to their real-space positions. The recon-space 
positions used are those which correspond to the number of iterations required to recover the known value of (3. The left panel illustrates the 
absolute displacement in km s“ ^, whereas the right panel shows the displacement of halos normalized by the difference between real-space and 
redshift-space positions. As such, a normalized displacement of -1 corresponds to a halo lying at its redshift-space position, and a normalized 
displacement of 0 corresponds to a halo returning to its real-space position. The solid blue curve corresponds to all halos, the turquoise dotted 
curve corresponds to halos with known velocities greater than 150 km s“ ^, and the dotted red curve corresponds to halos with known velocities 
less than 150 km s“ ^. Halos which have a normalized displacement less than -1 have moved away from their initial positions due to lying within 
triple-valued regions. Halos which have a normalized displacement greater than 0 have moved to their initial positions and beyond. 


their peculiar velocities). In this way the halos which lie in 
redshift-space have a normalized displacement of —1, and 
those which have returned to real-space have a normalized 
displacement of 0. With each successive iteration tracers 
would all ideally return to 0, i.e. return to the position they 
had in real-space. Halos which have a normalized displace¬ 
ment more negative than —1 have moved away from their 
initial positions. We argue below (§B) that these are halos 
in triple-valued regions. From this figure it can be seen that 
the typical tracer does in fact return to its real-space posi¬ 
tion. However, the distinction between low and high pecu¬ 


liar velocity tracers is more pronounced. The tracers with a 
low peculiar velocity have a symmetric distribution of recon¬ 
struction error. The halos which initially have a high veloc¬ 
ity, however, have a more skewed distribution which likely 
results in part from triple-valued regions. 

Having reconstructed the density field in recon-space, 
we return to the question of how well this field, once 
smoothed, can predict peculiar velocities. Fig. Al shows the 
results of this comparison. The fitted slope obtained from 
the reconstructed density field is systematically offset in 
comparison to those obtained from the real-space halo den- 
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sity field. The key result is that predicted velocities derived 
from the reconstructed smoothed halo fields yields an un¬ 
biased value of the slope, i.e. (3, when a Gaussian kernel 
of 4 h~^Mpc is used. This value is smaller than the value 
of ~5 h~^Mpc found for halos in real-space, because the 
recon-space density field is effectively pre-smoothed (in the 
redshift direction) by the ~ 300 km s“^ reconstruction error 
discussed above. 


APPENDIX B: ITERATIVE RECONSTRUCTION 
NEAR TRIPLE-VALUED REGIONS 

One problem with any reconstruction method is correcting 
for the fact that mapping from redshift to position is multi¬ 
valued, i.e. there can be several distances which map to the 
same redshift. 

To give an illustration as to how this affects the iterative 
reconstruction method, we have created a simple toy model 
in Fig. B1 which shows the velocity field around a spheri¬ 
cal overdensity. One can see that for any object with an ob¬ 
served redshift in the range ~ 100 to ~ 400 km s“^, there 
are multiple distances that correspond to that given redshift. 
As a concrete example, consider a galaxy at a redshift of 300 
km s“^. We will consider a simplified version of our recon¬ 
struction scheme in order to illustrate the general behaviour. 
In this simplified version, the value of 13 is held fixed, i.e. it is 
not increased adiabatically, so the amplitude of the peculiar 
velocity curve is fixed and the curve cz{r) is also fixed. 

In the first step of this iterative process, we place the 
galaxy at a distance corresponding to its observed redshift: 
300 km s“^ at position r(0) in Fig. Bl. At this location the 
predicted peculiar velocity is approximately —190 km s“^, 
so the object is moved to its first reconstructed position of 
490 km s“^ = 300 km s“^ — (—190 km s“^) at position 
r(l). At this position it has a predicted peculiar velocity 
of approximately —20 km s“^, and so this predicted pecu¬ 
liar velocity is again subtracted from its observed redshift 
of 300 km s“^, resulting in a reconstructed distance of 320 
km s“^indicated r(2). In this example, each successive loop 
brings the reconstructed position closer to the outer redshift 
solution of approximately 360 km s“^ indicated r((X)). The 
same behaviour occurs for redshifts between 250 and 400 
km s“^, the final position converges to the third (outer) of 
the three values. For redshifts between 100 and 250 km s“^, 
in contrast, the recon position converges to the first of the 
three values. One can see that the outer solutions are attrac¬ 
tors, and the middle solution is unstable. 

The above toy model is for a test particle (e.g. galaxy) 
moving in a fixed potential well (e.g. supercluster). It can 
also happen that two similar-mass galaxies (or two clusters) 
may each respond to the other’s gravity. In this situation, 
the distances of the two objects can “leapfrog”, i.e. can ex¬ 
change the order of their distances in reconstructed space. 
This leads to a sign change of the predicted peculiar velocity, 
which the leads to another “leapfrog” and a subsequent sign 
change in predicted peculiar velocity. Hence some objects 
show oscillatory behaviour in the reconstructed distance. 

Note however that the simple model described above 
and in Fig. Bl is an over-simplification of our reconstruction 
process. Of course, while in reality the density field used is 
a smoothed Gaussian random sphere not a spherical top-hat 



Real and Reconstructed distances (knn s ^ ) 

Figure Bl. Simplified reconstruction of one test object near a fixed 
potential well. In a Universe without structure, Hubble’s law (blue 
solid line y = x) would be sufficient to convert from redshift to 
real-space positions. In a Universe containing a spherical top-hat 
overdensity, however, the relation is complicated by infall towards 
the well (green solid curve). An object with an observed redshift in 
excess of 300 km s“ ^ can be seen to have originated from three pos¬ 
sible real-space positions (vertical red solid lines). The reconstruc¬ 
tion process starts by placing the particle at it’s observed redshift, 
predicting its velocity, and then subtracting that prediction from the 
original observation (the dotted lines parallel to the Hubble law line) 
to get a new predicted position. In most situations this new posi¬ 
tion will be beyond one of the possible real-space positions where 
the predicted infall velocities are small. Repeating the process here 
results in a much smaller velocity offset which can be projected 
(dot-dashed lines) back to the observed velocity to find the next re¬ 
constructed position. The iterative reconstruction converges in most 
situations to one of the two outer possible real-space positions. Re¬ 
constructed positions never converge, however, to the central possi¬ 
ble real-space position. 

overdensity, qualitatively the behaviour around overdensities 
is similar. More importantly, our iterative scheme is more 
complex than the one outlined above. First, we “adiabati¬ 
cally” increase (3* at each iteration. Consequently, the ampli¬ 
tude of the predicted peculiar velocity field increase slowly 
and so the size of triple-valued regions in redshifts space 
grows slowly. Second, we reduce oscillatory behaviour by 
averaging the positions of several successive iterations. Thus 
the convergence is better-behaved than the toy model de¬ 
scribed above. But qualitatively, the effect remains the same: 
as a result of preferentially placing objects towards the pe¬ 
riphery of triple-value regions, the density field is effectively 
“smoothed” more in the regions of highest density. 
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